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1 Introduction 

In the late eighties, Turbiner and Ushveridze [1] discovered some cases 
where a finite number of eigenvalues (and the associated eigenfunctions) of 
the spectral problem for the Schrodinger operator 

Hi> = Ei>, H = ~ + V(y) ,yeRorR + (1) 
dy 2 

can be found explicitely. The corresponding problems have been called quasi- 
exactly solvable (QES). Since that first step, QES equations have been classi- 
fied [2] according to their relation with the finite-dimensional representations 
of the Lie algebra sl2, R). Indeed QES Schrodinger equations as given in (1) 
can be written as 

d 2 d 

H(j> = E<j>, H = P4 (x)—+p 3 (x)—+p 2 (x) , xeRor R + (2) 

through ad-hoc changes of variables and functions [2] 

x = x(y) , ip = exp(x)<j), (3) 

if Pj(x)(j = 2,3,4) refer to polynomials of order j in x. The Hamiltonian H 
in (2) can also be expressed in terms of the first-order differential operators 

2 d 

]+ = -x — + 2jx, 
dx 

3o = x±-3, (4) 
d 1 

satisfying the si (2, R) commutation relations i.e. 

[jo,3±\ = ±j±, (5) 

\j + ,j.] = 2j , (6) 

the Casimir operator of this structure being C = j+j- + jl —jo- The operator 
H is then 

h= E O'^ + E c % (7) 

fi,v=±,0,fi>v /i=±,0 
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where the coefficients d 1 ", & 1 are arbitrary real numbers. 

The crucial point in order to relate the operator (7) to QES problems is 
the introduction of the nonnegative integer 2j in (4). Indeed, the generators 
of s/(2, R) as written in (4) preserve the space of polynomials of order 2j 



and so do the Hamiltonian (7). Searching for the eigenvalues of (1) is thus 
limited to the diagonalization of (7) in the (2j + l)-dimensional space P(2j). 
It is a straightforward problem leading to the knowledge of the eigenvalues 
Ek(k = 0, 1 , . . . , 2 j ) as well as the corresponding eigenfunctions ipk (cf. (3) 
where (fik belongs to P{2j)) and ensuring the quasi-exactly solvability of the 
original equation. 

We propose in this paper to take a new look at this problem through the 
consideration of the so-called polynomial deformations [3] of sl(2, R) i.e. of 
the structures characterized by the following commutation relations 



where p n (Jo) stands for a polynomial of order n in the operator J . 

More precisely, in Section 2, we show how to introduce in a natural man- 
ner the polynomial deformations (10) inside the operator (7) (n will then be 
limited to 3). In Section 3, we study the finite-dimensional representations 
of these polynomial deformations while Section 4 is devoted to their finite- 
dimensional differential realizations. In Section 5, we analyze two examples 
i.e. the sextic oscillator (Subsection 5.1) and the second harmonic genera- 
tion (SHG) problem (Subsection 5.2). Finally, we give some conclusions in 
Section 6. 



2 Polynomial deformations of s/(2, R) inside 
QES problems 



Instead of considering the operators (4) as expressed in (7), let us in- 
troduce the following operators (so defined for natural reasons in connection 
with their respective raising, diagonal or lowering characteristics) 



P(2j) = {l,x,x 




(8) 



[<^0> — =t^±, 

[J+,J-\ =Pn(Jo), 



(9) 

(10) 



J+ = + c +0 3+3o + C + J+, 



(11) 
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J = c+- 3+ j„ + c 00 f +c j , (12) 

J_ = c°-j j_ + c—ji + c-j_, (13) 
so that H simply writes 

H = J+ + J + J- (14) 

As we will see, restoring the linearity inside (7) such that it becomes the 
combination (14) will have the consequence of introducing nonlinearity inside 
(6) so that we will be concerned with the algebra (9)-(10). Indeed asking for 
the relations (9) to be satisfied with the operators (11)-(13) and the relations 
(5)- (6) lead to two cases only i.e. either 

C++ = c— = 0, c+- = c 00 , c° + c 00 = 1 (15) 

or 

c ++ ^ o, c— ^ 0, c+° = c+ = c°- = c" = 0, c+- = c 00 , c° + c 00 = i. (16) 

The relation (10) is then ensured with the respective polynomials p n= 3(Jo) 

p 3 (J ) = 4c+°c°- J 3 + 3(c°-c+ + c+°c- - c+°c°-) J 2 

+ [2c+c- - c+°c- - c°-c+ + c+°c°-(l - 2j(j + 1))] Jo (17) 

+j(j + l)(c +0 c°- -c+V -c°-c+) 

or 

p 3 ( Jo) = 8c++c-((2j 2 + 2j - 1) J - 8J 3 ). (18) 

Notice that in these expressions, c+~(= c 00 ) has been put equal to zero 
without loosing generality (cf. the Casimir operator of sl(2,R)). Moreover, 
the relation (18) refers to the Higgs algebra already intensively visited [3]. 

We can thus consider the QES Hamiltonians as linear combinations of 
operators generating the following polynomial deformation of sl(2,R) 

[J ,J ± ] = ±J±, (19) 

[J+, J_] = a J 3 + (3Jl + 7 J + 5, a,P,~f,8eR (20) 

taking account of the two possibilities (17) and (18). Notice that the Casimir 
operator of this deformed algebra is 

= J + J - + 4 J o + ( 3 - 2 } J ° + ( 4 " 2 + 2 } J ° + ( 6 " 2 + 5) J ° 
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The next step will be the determination of the finite-dimensional repre- 
sentations of the algebra (19)-(20) denoted in the following by sl^(2, R), the 
upper index referring to the highest power of the diagonal operator. 



3 Finite-dimensional representations of s/( 3 )(2, R) 

As stated in the Introduction and in relation with the possible diagonaliza- 
tion of if, we are interested in the finite-dimensional (= 2J+1, J — 0, |, 1, ...) 
representations of sl^(2,R), only. We thus consider kets of type \J,M> 
with M running from —J to J and such that 

M 

J | J,M >= (— + c) | J,M >, (21) 
q 

J+\ J,M >= f(M) | J,M + g>, (22) 

J_ | J,M >=g(M) | J,M-q >, (23) 

where q is a positive integer and c a real number. The relations (21)-(23) are 
such that (19) is satisfied. In order to ensure (20), we have to impose 

f(M - q)g(M) - f(M)g(M + q) = a(j + c) 3 + /3(j + cf 
M 

+ c) + 5, M — -J, J. (24) 

Moreover, we have to take account of the dimension of the representations, 
leading to the constraints 

f(J) = /( J - 1) = ... = /( J - q + 1) = (25) 

and 

g(-J) = g{-J + 1) = ... = g(—J + q-l)=0. (26) 

Being interested in the highest weight representations, we obtain from (24) 
and (25) the following result 

f(J — (k+ l)q - l)g{J -kq-l) = (k+ !){«(— - + c) 3 + - + c) 2 
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— - + c)+5- l[3a(— - + c) 2 + 2/3(— - + c) + 7 ]fc 
q 2 q q 

+ l [3a (^r + c) + ® k ^ 2k ^ k2{ - k + 1)} 



(27) 



where / = 0, 1, q — 1 and k = 0, 1, 



2J-d-l 



. The nonnegative integer d 



introduced in the last formula has to take specific values according to / but 
also to J. These values are summarized in the following table, n denoting a 
nonnegative integer. 



Table 


/ = 


Z = 1 


I = 2 




Z = q- 1 


J=(qn)/2 


d = 


d = q — 1 


d = q — 


2 




d = 1 


J=(qn + l)/2 


d = 1 


d = 


d = q — 


1 




d = 2 


J =(qn + 2)/2 


d = 2 


d= 1 


d = 






d = 3 


J= (qn + q- l)/2 


d = g-l 


d = q-2 


d = q — 


3 




d = 



Taking care of these values, we also have to constrain the real c in (21) and 
(27) through the conditions (26). This leads to equal to zero the expression 
inside the brackets { . } in (27) for k = 2J ~ d ~ / or, in other words, to consider 
the q equations on c 



1 



1 



a[c 6 + — (d - l)c 2 + —(J 2 - J(d + l)+l 2 -dl + d 2 )c +—(2J-d- l)c 



2q 



2q 



1 



+ ^(2J 2 (d - I) - 2J(d 2 - I 2 ) + d 3 - d 2 l + dl 2 - I 3 ) + ^(l 2 -d 2 + 2J(d 
+[3[c 2 + -(d-l)c+ -^(J 2 - J(d + l) + d 2 -dl + I 2 ) + — (2J -d-l)} 



q 
l 



3q 



6g' 



+ 7 (c + -(d - Z)) + 5 = , Z = 0, 1, q - 1. 
iq 



Excluding the trivial case (a, (3, 7, 5) = (0, 0, 0, 0), we notice that these equa- 
tions reduce to 



ac(c 2 + J{J + 1)) + [3(c 2 + -J{J + 1)) + 7 c + 5 = 0, 



for q — 1 while for g = 2, we are led to either 

a = 0^/5 = 0,c = 

6 
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(29) 



(30) 



or 

^°^ 5= 3^-27^' C= -^ (31) 



if J = n and to 



a(3c 2 + -V(J + l)-i) + 2/3c + 7 = 0, (32) 
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a(c 3 + i J( J + l)c) + /3(c 2 + J( J + 1)) + 7 c + 5 = 0, (33) 

if J = n + \. Two other values of q are also available namely q = 3 and 
g = 4. We respectively obtain 

^°' 7= 3^-9 aJ + 9 a ' 5= 3^-27^' C= -^ (34) 



if J=f, 

1 '-v = .... _ 

3a 9 9 9 ' 3a 27a 2 ' 3a 



/3 2 1 , 2 2 1 /3 7 2/3 3 



iiJ= 3n±l 

/5 2 1 .2 1 ,^1 A 2/3 3 

a =t U ,7 = «J aJ H — a,o = n,c = (do) 

T ' ' 3a 9 9 9 ' 3a 27a 2 ' 3a V ; 

if J = 3i t± 2 -, these three contexts being associated with q = 3. We also have 

, n /5 2 1 o 3 /? 7 2/5 3 /3 

U,7= aJ H a,o = ~,c= (37) 

T ' ' 3a 16 16 ' 3a 27a 2 ' 3a V ; 

if J = 2n and 

^ 1 i* 1 ,1 x Pi 2/3 3 (3 
a ^°^ = ^--6 aJ -8 aJ+ 8 a ' 6 = ^-^ C= -M (38) 

if J = 2n+l, these two cases being related with q — 4. The two other systems 
related to g = 4 i.e. those corresponding to J = 2n + \ and J = 2n + | are 
incompatible ones as are those related to q > 4. 

Let us also notice that some of these representations are reducible. For 
instance, if we consider the case of the usual sl(2,R) algebra, corresponding 
to a = 0, (3 = 0, 7 = 2, 5 = 0, we obtain 

q = 1 -> c = (39) 
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and 

q = 2->c = 0,J = n. (40) 
The first case (39) is associated to (see (27)) 

f(J-k- l)g(J — k) — (k + 1){2J - k} , k — 0, 1, 2J (41) 

or, in other words, to 

f(M - l)g(M) = (J - M + 1)( J + M), M = - J, J. (42) 

We recognize in (42) the well known result of the angular momentum theory 
[4] subtended by this sl(2,R) algebra. The second case (40) corresponds to 

f(J -2k- 2)g(J - 2k) = (k + 1){ J - k} , k = 0, 1, J, (43) 

f(J - 2k - 3)g( J - 2k - 1) = (k + 1){ J - k - 1} , k = 0, 1, J - 1, (44) 
or, in other words, to 

/(M - 2)#(M) = i( J - M + 2)( J + M), M = -J, - J + 2, J, (45) 

f{M-2)g(M) = I(J-M+1)(J+M-1), M = - J+l, - J+3, J-l. (46) 

It is then clear that the representation ( J = n, q = 2) is in fact the direct sum 
of the two (irreducible) representations (J = ^^^q = 1) and (J = = 1) 
(the eigenvalues of the Casimir being equal). 

4 Finite-dimensional differential realizations 

of sl^(2,R) 

We now turn to the construction of the differential realizations (expressed 
in terms of the real variable x) of the algebra (19)-(20). In correspondence 
with (21)-(23), the generators of sl^ z \2,R) have the following forms 

J+ = a?F(D), (47) 
J =l(D-J) + c, (48) 
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J. . G(D)£. (49) 

and the basis {| J, M >, M — — J, J} stands for the space P(2J) (cf. (8)) 
of monomials {x J+M , M — —J, J}. Moreover, D is the dilation operator 

D = x A. (50) 



By remembering that 



and 

d q D! 



a; 9 



the relation (20) gives 



F(D - q)G(D - g)^ - F(D)G(D)^_±iL = 
«[-(£>- J) + c] 3 + /?[-(£- J) + c] 2 

+7[^(L>- J) + c] + 5. (53) 

Let us discuss this condition within the (a ^ 0)-case first and the (a = 0)- 
case second, a) In order to avoid singularities, we thus impose, when a^O, 

F(D)G(D) = -^ip^yyP + Ai)(D + A 2 )(D + h)(D + A 4 ) (54) 

and obtain the following system on the real unknows Ai, A4 

4 3 

Ai + A 2 + A 3 + A 4 = 4gc + 2g-4J + --g, (55) 

6 a 



AiA 2 + + A1A4 + A 2 A 3 + A 2 A 4 + A3A4 = Qq 2 c 2 + 6q 2 c - 12gJc 

+A-q 2 c + 6J 2 - 6gJ - 4-gJ + q 2 + 2-q 2 + 2-q 2 , (56) 
a a a a 
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A1A2A3 + A1A2A4 + A1A3A4 + A2A3A4 — 

4g 3 c 3 - 12q 2 Jc 2 + 6q 3 c 2 + A-q 3 c 2 + 12qJ 2 c - Y2q 2 Jc 

a 

+2q 3 c - 8-q 2 Jc + A-q 3 c + A-q 3 c - 4J 3 + QqJ 2 - 2q 2 J 

a a a 

+4 d qJ 2 _ 4—q 2 J + M g 3 _ A l q 2 j + 2 7 g3 + 4 i g3> (57) 
a a 3a a a a 

Let us recall that we are interested in finite-dimensional (=2 J+l) realizations 
only. This means that the cases g = 1 and q = 2 are the only possibilities in 
accordance with 

q = 1 _> Ai = 1,A 2 = -2J (58) 

and 

g = 2 -> Ai = 1, A 2 = 2, A 3 = -2 J, A 4 = —23 + 1. (59) 

In the first context (q — 1), the equations (55) and (56) fix A 3 and A 4 as 
follows 

^3 = 2c+i-J+^ + i \ 4 = 2c+ 1 --J+^- e - (60) 
2 6a 2 2 6a 2 

with 

f ^ l+ l ^i.-U(J + l)-^- l 4 l ic-sl (61) 
9a 2 3 a a 

while the equation (57) coincides with (29). In the second context (q = 2), 
we are led to (31) supplemented by 

B 2 a t9 a T ct 

b) The case a = is more simple and has already been analyzed [5]. For 
self-consistency, we recall the main results i.e. 

F(D)G(D) = - J^^L^D + A0(D + A 2 )(D + A 3 ) (63) 

where the only possible finite-dimensional (=2J+ 1) realization is associated 
with q — 1 and corresponds to 

A x = 1, A 2 = -2 J, A 3 = - J + 3c + ^ + 1 (64) 
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the real c being fixed through 



(3c 2 + 1C + 5+^J(J + l) = 0. (65) 

Now that both cases have been considered, let us conclude this Section 
by noticing that some realizations (namely the ones corresponding to q = 2 
without the condition (62) and the ones associated with q — 3, 4) are missing 
with respect to the representations developed in the previous Section. Indeed, 
the relation (53) we have imposed is more constraining because it is a relation 
between operators independently of the basis on which they are supposed to 
act. If we take account of this basis i.e. in this case P{2 J) = {x J+M , M = 
— J, J}, we can recover all the cases previously discussed. For example, in 
the context q — 3, J = |, we can consider 

J+ = -\!(-\)^(D - l)(D - 2)[D - 3), (66) 

J - - < 68 > 

with 

f(-\)9i\) = f ■ (69) 

It is then easy to convince ourselves that these operators generate sl^(2, R) 
with 7 = ^ — ^ and 5 = — but on the space -P(3) only (the relation 
(53) being trivially not satisfied except on this space). However it has to 
be stressed that the relations corresponding to (66)-(68) but in the general 
context become really heavy when the value of J increases. 



5 Two examples 

We first consider the prototype of QES systems i.e. the so-called sextic 
oscillator [1,2] and then turn to a more physical example: the SHG problem. 
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5.1 The sextic oscillator 

This system is characterized by the following potential 

V(y) = a 2 y 6 + 2aby A + (6 2 - 2ap - 8aj - 3a)y 2 (70) 

with 0), b e R and p = 0, 1 while j is the quantum number appearing in 
(4). With 

x = y 2 (71) 

and 

* = -/<r + |-i>< < 72 > 

we can be convinced of its QES character via the form (7) (up to a translation) 
H = J + + J + J- + 2fy» + 46j + 6 (73) 

and 

c°~ = -4, c + = -4a, c° = 46, c" = -(4j + 2 + 4p), (74) 

the other c's being equal to zero. Without loss of generality, we can put 
b — | (in order to recover (19)) and obtain, through (11)-(13), the relation 
(20) with 

a = 0,P = 48a, 7 = 32a(p + j),6= -16a j(j + 1). (75) 

The case q — 1 is thus the only one to be available. We actually have 

Jo I J,M >= (M + c) I J,M >, (76) 

J+ I J,M>=f(M) I J,M + 1 >, (77) 
J_ I J, M >= g(M) I J, M - 1 >, (78) 

with 

f{M-l)g{M) = (J-M + l)(J + M)(A8ac+lQaM+lQap + lQaj-8a). (79) 
Moreover, the parameter c is fixed according to 

c = ~\ip + J) ± \\[iP + J) 2 ~ 3 J( J + 1) + 3j(j + 1) (80) 
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leading to constrain J through 



J <-- + -^Qj(j + l) + 12(p + j)i + 9 (81) 

in order to ensure the reality of c. Because the space P(2J) is preserved, we 
just have to equal to zero the determinant of the following matrix 

2J+1 o 

M= £(E + J-k-c-^-j + -)e k , k 
k=i z 4 

2J 2J 

" E 9(~ J + k)e k)k+1 - f(~J + k- l)e k+1>k (82) 

k=l k=l 

in order to find the energies. In the matrix (82), the notation e k ,i stands for 
a (2 J + l)-dimensional matrix where 1 is at the intersection of the k th row 
and the I th column, the other elements being 0. For example, when j — |, 
we have 

J = 0,i (83) 
according to (81) while the relation (80) gives 



c — i-i ± iv & ' + i ) "--i (84) 

if J = and 

c = or c = --(2p+ 1) (85) 
o 

if J = |. In the case (84), the energies are 

E = c + - -\ > E = 0.0428932; 0.0562872; 1.1103796; 1.4571067 (86) 

2 4 

and in the case (85), the resolution of the vanishing determinant associated 
with (82) leads to 

E = -± - Vl + 32a; E = -± + 96a (87) 
4 2 4 2 

if c = and 

E = — ± -VI - 32a; E = -± - 96a (88) 
12 2 '42 y J 
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if c = — |; c— — 1. Only the values given in (87) correspond to the previously 
obtained ones [2] . This is indeed a general result that the known energies [2] 
are recovered through our approach when c = 0, J = j. The other contexts 
(J < j) or (J — j,c — — |(p + j)) lead to supplementary new values of the 
energy. 

Let us analyze more deeply this result by going to the differential realiza- 
tion (47)-(49) i.e. 

J+ = xF(D), (89) 
J = D - J + c, (90) 

J - - < 91 > 

where (cf. (63) and (64)) 

F(D)G{D) = -16a(D-2J)(D - J + 3c + - + p + j). (92) 

In order to preserve the space P(2J), let us make the choice (without loss of 
generality, this freedom being due to the fact that s/( 3 )(2, R) is defined up to 
an automorphism 



G{D) = -4(£>- J + 3c+^+p + j). (93) 
In that case, the Hamiltonian (73) is realized as 

H = -4x-^r + [4ax 2 + x -4(-J+ 3c + ^ +p + j)]-^- 
dx z 2 ax 

-8a Jx + | + 4+ J - J + c - ( 94 ) 

This form is analog to (2) and we propose to write it in the Schrodinger form 
(1) through the changes (3) i.e. 

x = y 2 , (95) 

ij = exp(J {^-2ax+^(J-3c-p-j)^)dx)<j>. (96) 
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The potential obtained in this manner is given by 

V(y) = a 2 y 6 + ^ a y 4 + (y^ — 6aJ — 2aj — 6ac — 2ap — 3a)y 2 

+l(j-J-c) + (J-3c-p-j)(J-3c-p-j + l)\. (97) 
2 y l 

Compared with (70), this expression actually reduces to the sextic oscillator 
potential iff c = and J = j. For other values of the parameters (i.e. the 
already cited (J < j) and (J = j, c = — |(p + j)), the new eigenvalues of 
the problem (appearing for j — \ in (86) and (88)) do correspond to another 
model, namely the radial sextic oscillator as shown by (97). 

5.2 The second harmonic generation 

This nonlinear optical process as well as others such as coherent spon- 
taneous emission and down conversion [6] can be described by the following 
effective Hamiltonian 

H = a\a x + 2a\a 2 + g{a\a\ + (a\) 2 a 2 ) (98) 

with cubic terms in the (independent) bosonic creation and annihilation op- 
erators. The Hamiltonian (98) has already been recognized [7] as a QES 
model. We confirm such a result by using the technique developed in Section 
2. Indeed following Section 2, we propose to define 

J+ = a\a\, (99) 

Jo = ^(4 Q 2 - a\ai), (100) 

J- = (a\) 2 a 2 (101) 
such that the algebra (19)-(20) is satisfied with 

a = 0, (3 = —12, 7 = 0, 5 = 7-iV 2 + N. (102) 

o 

In the last expression, iV is the invariant 

N = a\ai + 2a\a 2 (103) 
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satisfying 

[N, Jo] = [N, J±] = 0. (104) 

Once again, the values (102) are typical of the q = 1-representation only so 
that the relations (76)- (78) are the ones to be taken care of but with 

f(M - l)g(M) = (J - M + 1)( J + M)(-12c - 4M + 2). (105) 

The parameter c is fixed according to 

c2 = 4 J(J + 1) + i Ar2 + B /v < 106) 

and its reality is ensured if 



J<-\ + \f-Ni + N + l. (107) 

It is then possible to determine the energies by putting to zero the determi- 
nant of a (2 J + 1) by (2 J + 1) matrix analog to (82) as well as it is possible 
to determine them through the differential realization (89)-(91). In this case, 
we have 

F(D)G(D) = A(D-2J)(D- J + 3c+i) (108) 

and choosing 

G(D) = 1 (109) 

the Hamiltonian (98) becomes 

H = N + Agx z ^ + g{l + 12(- J + c + \)x 2 )4~ 
ax z 2 ax 

+4g(-J + 2J 2 -6Jc)x. (110) 
With the respective changes of variables and wavefunctions 

x = -—. - v i> = exp{- f (_L + -( C - J)-)dx)<j> (111) 
gy 2 J 8x A 2 x 

we can put (110) on the Schrodingerlike form (1) with 

V(y) = (-y« + \g\c -J- \)y 2 + ( J + 3c) (J + 3c + 1)1 + N. (112) 
lb 2 2 ?/ z 
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Once again this is typical of a radial sextic oscillator and the QES character- 
istics of the second harmonic generation is thus proved. The determination 
of the energies is then possible without any difficulty [2]. For example, when 
N = 4, we have 

J = 0^£ = 4, (113) 
J = - -> E = 4±#v / 2v / l9, (114) 

J = 1 _> E = 4,4±4:g. (115) 

Only the values (115) correspond to known ones, the values (113), (114) 
coming from other models. This is a general result in the sense that the 
energies of SHG are the ones of the Schrodinger potential (112) with 

N 

c=J-j (H6) 

and 

J-$-J-*i± (117) 
according to even or odd values of N. The SHG potential thus writes 

V(y) = ^y 6 - 9 j(2N + 3)y 2 + N (118) 

and exactly coincides with the potential (42) of Ref. [7]. In terms of the 
operators (4) (with j = J), this gives 

1 N 

H = 9u(jo + ^ - j) - ±93- + N (119) 



if N is even and 



if N is odd. 



H = gj + {j ~ \ ~ j) - 4gj- + N (120) 
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6 Conclusions 



We have developed a general method based on the polynomial deforma- 
tions of the Lie algebra sl(2, R) in order to exhibit the QES characteristics 
of a Hamiltonian. We have applied this method to two examples: one more 
theoretical -the sextic oscillator- and one more physical -the second harmonic 
generation-. In both finite number of energies as well as eigenfunc- 

tions are determined through the finite-dimensional representations -or, in 
an equivalent way, through the realizations- of these polynomial deforma- 
tions. Some of these energies (and eigenfunctions) were previously known, 
not the others. It seems, through the analysis of the two examples, that 
these previously unknown energies do not correspond to the same model but 
to another one being closely related to the first one. In some cases, the com- 
parison between these new and old models could be interesting, with respect 
to experimental data in particular. 

The main advantage of the method we have proposed is that it can be 
systematically applied to any Hamiltonian written in terms of a raising and 
a lowering operator. Numerous physical Hamiltonians are of that type. In 
particular, we plan to analyze one of them: the so-called Lipkin-Meshkov- 
Glick Hamiltonian [8] of specific interest in nuclear physics. Being based on a 
(a 7^ 0) polynomial deformation, its analysis is more delicate but also richer 
in the number of available representations [9]. 

The main drawback is that it is limited to sl(2,R) and its deformations 
when we know that some QES models need more extended algebras. How- 
ever, the method we have presented here can also be generalized to these 
extended algebras as well as superalgebras. We also plan to come back on 
these points in the near future. 
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